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Abstract 



The low energy radiative capture process np — ► d^y provides a sensitive probe 
of the two-nucleon system. The cross section for this process is dominated by 
the isovector Ml amplitude for capture from the channel via the isovector 
magnetic moment of the nucleon. In this work we use effective field theory 
to compute the isoscalar Ml and isoscalar E2 amplitudes that are strongly 
suppressed for cold neutron capture. The actual value of the isoscalar E2 
amplitude is expected to be within ~ 15% of the value computed in this 
work. In contrast, due to the vanishing contribution of the one-body operator 
at leading order and next-to-leading order, the isoscalar Ml amplitude is 
estimated to have a large uncertainty. We discuss in detail the deuteron 
quadrupole form factor and 3 S\ — s Di mixing. 
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The cross section for radiative capture np — > d'-f of thermal neutrons has an important 
place in nuclear physics as it provides a clear demonstration of strong interaction physics 
that is not constrained by nucleon-nucleon scattering phase shift data alone. Effective range 
theory |],|2| uniquely describes the scattering of low-energy nucleons, yet fails to reproduce 
the measured cross section of <r expt = 334.2 ±0.5 mb (measured at an incident neutron speed 
of |v| = 2200 m/s) for np — > dj at the 10% level. In the effective field theory appropriate 
for very low momentum interactions || (i.e. without pions), this discrepancy is understood 
to arise from the omission of a four-nucleon-one-magnetic-photon operator that enters at the 
same order as effective range contributions. Conventionally, this discrepancy is attributed to 
pion-exchange-currents . The cross section for np — > d'-f at very low energies is dominated 
by the capture of nucleons in the ^So state, via the nucleon isovector magnetic moment, the 
amplitude for which we denote by Mly. This particular amplitude is much larger than 
other amplitudes for several reasons. First, initial state interactions give a contribution 
proportional to the large scattering length in the ^So channel, a*- = —23.714 ± 0.013 fm. 
Second, the Mly amplitude is proportional to the nucleon isovector magnetic moment, K\, 
which is much larger than the nucleon isoscalar magnetic moment, kq, which dictates the size 
of the one-body contribution to the isoscalar magnetic amplitude, Mis- Third, the capture 
from the 3 5*i channel that does proceed via the nucleon isoscalar magnetic interaction (the 
one-body contribution) must vanish at zero-momentum transfer as it is the matrix element 
of the spin operator between orthogonal eigenstates. Finally, the electric amplitudes, Ely 
for capture from the P-wave and E2s for capture from the 3 Si channel, are suppressed by 
additional powers of nucleon momentum or photon energy compared to Mly. 

While the Mis, E2s and Ely amplitudes are much smaller than Mly, measurements of 
spin-dependent observables can determine specific combinations of these amplitudes. Two 
such observables are the circular polarization of photons emitted in the capture of polarized 
neutrons by unpolarized protons, and the angular distribution of photons emitted in the 
capture of polarized neutrons by polarized protons. The circular polarization of photons 
emitted in the forward direction in the capture of polarized neutrons on unpolarized protons 
has been measured to be p^xpt = —(1.5 ± 0.3) x 10~ 3 . This value is consistent with 
previous theoretical estimates |J . An experiment that will measure the angular distribution 
of photons emitted in the capture of polarized neutrons on polarized protons is to be carried 
out at the ILL reactor facility [[| and results should be available in the near future. 

In this work, we calculate the Mis and E2s isoscalar amplitudes that contribute to 
np — > d'-f using the effective field theory (EFT) of nucleon-nucleon interactions without pi- 



ons, EFT(ft) } as detailed in ||, using KSW power counting fll~9| , |21~l| . A significant amount of 
progress has been made in the application of EFT to the two- and three-nucleon systems 
fi~0|j - f3~9fl during the past few years. A test of this formalism will be a comparison between 
these predictions for the strongly suppressed amplitudes in np — > d-f and the measured 
experimental asymmetries which constrain them. Calculations of these suppressed ampli- 
tudes using an alternative power counting are being performed by Park, Kubodera, Min and 
Rho [50]. Our work results from a challenge issued by M. Rho for the community to make 



predictions for these amplitudes ||0 



The amplitude for low-energy np — > d'-f is 
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where we have shown only the lowest partial waves, corresponding to electric dipole capture 
of nucleons in a P-wave with amplitude Xei v , isovector magnetic capture of nucleons in 
the ^So channel with amplitude Xmi v , isoscalar magnetic capture of nucleons in the 3 S\ 
channel with amplitude Xmi s , and isoscalar electric quadrupole capture of nucleons in the 
3 5i channel with amplitude Xe2 s ■ As we dimensionally regulate the divergences that appear 
in the effective field theory we keep explicit space-time dependence in the amplitudes shown 
in eq. ([1|), with n the number of space-time dimensions. U n is the neutron two-component 
spinor and U p is the proton two-component spinor. p is half the neutron momentum in 
the proton rest frame, while k is the photon momentum. The photon polarization vector is 
e( 7 ), and £(d) is the deuteron polarization vector. For convenience, we define dimensionless 
variables X, by 
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where 7 = y/M^B ~ 45.6 MeV is the deuteron binding momentum, with B the deuteron 
binding energy. 

By measuring certain observables of the np — > process the four amplitudes X E i v , 
X Mlv , Xmi s , and X E2s , can be determined or constrained. The simplest quantity to mea- 
sure is the total cross section for the capture of unpolarized cold neutrons with speed |v| 
by unpolarized protons at rest (the neutron velocity v is related to the momentum p by 
p = |Mtvv, where relativistic corrections have been neglected). In terms of the amplitudes 
given in eq. (fj) and eq. (fj) the unpolarized cross section is 
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where a is the fine-structure constant. The cross section for the capture of cold neutrons is 
dominated by Xmi v by several orders of magnitude and therefore a measurement of a does 
not constrain the other three amplitudes. 

A spin-polarized neutron beam incident upon a spin-polarized proton target enables 
spin-dependent observables to be measured, even without measuring the polarization of the 
out-going photon or deuteron. If the protons have polarization r] p and the neutrons have 
polarization t] n , along the direction of the incident neutron momentum, the spin-dependent 
capture cross section is 
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where 9 is the angle between the polarization axis and the direction of the emitted photon. 
Spin-averaging the expression given in eq. (||) over the initial nucleon spin states, (r/ n , r/ p ) = 
(±1, ±1) and integrating over all angles reproduces the spin independent cross section shown 
in eq. (|3]). From this one can define the angular asymmetry, 
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For systems with high polarization, measurement of this angular asymmetry constrains the 
small amplitudes. In the expressions for S^ 1 ' and that appear in eq. (f| we have neglected 
the small Mis, E2s and Ely amplitudes in the denominators. 

If the polarization of the out-going photon can be measured, then other spin- dependent 
observables can be considered. For a polarized neutron incident upon an unpolarized proton 
target, there is a different cross section for production of right-handed versus left-handed 
circularly polarized photons. Defining the asymmetry (9) to be the ratio of the difference 
to the sum of these cross sections, 



AU9)=r} n (P 7 (M1) + P 7 (£2)) cosfl + P y {El) sin 2 9 



(7) 



where 



P 7 (M1) 



\^2Re[XMi v Xl 



\X 



P 7 (£2) 



^/2Re[X M i v X E2s ] 



Mix 



\X 



Rg[Xmi v X e1i 
\Xmi v I 2 



(8) 



where we have again neglected the small MI5, Elg and Ely amplitudes in the denominators. 

The four amplitudes X E i v , X Mlv , Xmi 8 , and X E2s , can be computed with EFT(f). 
Power counting the leading order (LO) versus next-to-leading order (NLO) for a given am- 
plitude is straightforward and follows the well known power counting rules |4],|19|,^1]] . How- 
ever, power counting amplitudes relative to each other is not so straightforward. The reason 
for this is that there are two different kinematic scales for the capture of cold or thermal 
neutrons - the photon energy and the momentum of the incident neutron. While the ve- 
locity of the incident neutron is always assumed to be small, its finite value gives rise to an 
Ely amplitude, which for |v| = 2200 m/s is comparable to the subleading MI5 and E2$ 
amplitudes. 
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It is convenient to express the X amplitudes as a series in powers of Q; X = X(~V + 
+ X^ + • • • where Q ~ j/m^ is the small expansion parameter in the theory and 
superscripts denote the order in Q. The isovector Ml amplitude X M i v has been computed 
with EFT previously p|,p6| up to NLO. The amplitude starts at Q° in the power counting, 



m v = * (i - « (15o) t) , o) 

where )/2 is the isovector nucleon magnetic moment in nuclear magnetons, with 

k p = 2.79285, n n = -1.91304. While naively, X$ lv is of order Q°, numerically xf} ly ~ 20 
due to the large numerical values of both k± and a^ 50 ^. 

At order Q 1 there are contributions to Xmi v from insertions of the effective range pa- 
rameter and also contributions from a four-nucleon-one-magnetic operator, described by the 
Lagrange density 

C = e% (N T PiNf (N T P 3 N)Bi + h.c. , (10) 

where B = V x A is the magnetic field operator. Pj and Pj are the 3 Si and 1 Sq spin-isospin 
projection operators respectively, with 

Pi = -77T°2(Ti r 2 , Pi = -7=02 r 2 Ti . (11) 



The NLO contribution to the amplitude is found to be [[4"l,f26| 
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where r = 2.73 ± 0.03 fm is the effective range in the ^So channel and pd = 1.764 fm is 
effective range in the 3 Si channel, p is the renormalization scale, and the //-dependence 
of ^L]_ yields a renormalization scale independent amplitude, by construction [|4|,p6"|. For 
convenience we choose p = m^. As Xmi v is the dominant amplitude for the capture process, 
= 7.24 fm 4 from the unpolarized cross section [|]] in eq. (|3|). 

The cross section for any finite incident nucleon momentum has a contribution from 
isovector El capture. Recently, a N 3 LO calculation of this amplitude has been performed 
plfl for non-zero energy capture. At N 3 LO there are contributions from the effective range 
parameter and from P-wave initial-state interactions which are found to be small. Neglecting 
the P-wave initial-state interactions, the amplitude is found to be, up to N 3 LO 

x E1 „ = -!p!^(i + i TO + 5 7 v3 + l 7 3p3) . (is) 

Capture from the P-wave introduces the factor of the external nucleon momentum, |p|, 
forcing the amplitude to vanish at threshold. The powers of jpd that appears in the ampli- 
tude are consistent with the deuteron S-wave normalization factor 1 / ^1 — ^pd that arises 
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in effective range theory. For moderate incident momenta, where |p| ~ Q, the LO Ely 
amplitude is of order Q~ x , and dominates the isovector Mly amplitude, which starts at 
Q°. However, for smaller incident momentum, the Ely amplitude becomes less important. 
If we take |p| ~ Q 2 , the Ely and Mly amplitudes are of the same order in the counting, 
however, for the neutron incident velocity of 2200 m/s, numerically |p| ~ Q . 

In the zero recoil limit, the matrix element of the nucleon magnetic moment operator 
between the deuteron and nucleons in the 3 S\ channel, contributing to Mis, is the matrix 
element of the spin operator between orthogonal eigenstates states of the strong interaction 
and thus vanishes. This leads to = at LO (Q°) and further, the contribution from 

the one-body operator at NLO (Q l ) also vanishes. However, at NLO there is a contribution 



from a four-nucleon-one-photon two-body operator defined by the Lagrange density [[J.^l 



C = -e *L 2 ie ijk (iV T iW) f (iV T P, iv) B fc + h.c. . (14) 



At NLO the deuteron magnetic moment is found to be |21] 



= 2Mjv ( S + Kn+ U n ) ■ (15) 

Reproducing the experimentally observed value of the deuteron magnetic moment requires 
that, at this order ^L 2 (m w ) = -0.149 fm 4 , which is significantly smaller than the 

naively estimated size of ~ 1 fm 4 . This two-body interaction contributes to X^ ls , and at 
NLO we find 

*fil = V2^L 2 (/,- 7 ) 2 . (16) 

While formally the leading contribution, the smallness of ^L 2 suggests that the contribution 
given in eq. (|T6| ) might not dominate over higher order terms, and the Ml $ amplitude might 
not be predicted well by EFT(^) at this order. 

To make this more concrete, one can imagine a higher dimension four-nucleon-one-photon 
local operator that gives rise to a contribution of the form 

AX Mls ~ V2^- *L X (n - 7 ) 3 Up 2 + 7 2 ) + (p' 2 + 7 2 )] , (17) 

between states with nucleon momentum p and p', since ^Lx ~ (// — 7) -3 - This object makes 
a vanishing contribution to the magnetic moment of the deuteron, while making a non-zero 
contribution to the rate for capture from the 3 Si channel 

AX^-^^^-T) 3 • (18) 
Ztt 

One naively expects ^Lx ~ 1 fm 6 , which would make such a contribution approximately 60% 
of the amplitude in eq. fliCf ). This relatively large uncertainty in the MI5 matrix element is 
consistent with previous calculations of this quantity M, and particularly the most recent 



(preliminary) work of Park, Kubodera, Min and Rho |40], where they find that different 
treatments of the short-range component of the interaction leads to an approximate 60% 
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uncertainty. At higher orders, there is a contribution from the one-body operator due to 
the finite energy release of the capture process. Naively, this contribution is much smaller 
than the expected contribution from higher dimension operators, as estimated in eq. ([TED, 
and so we do not consider it further. 

The E2s amplitude is dominated by local operators that convert 3 Si states to 3 D\ states 
and vice versa. However, the relatively slow convergence in this channel requires that the 
calculation be performed to higher orders so that a meaningful estimate of uncertainties is 
possible. In previous works [|J the deuteron quadrupole form factor and the 3 Si — 3 Di mixing 
parameter e\ P2"| , |43| were computed up to NLO. Presently, we compute E\ up to N 3 LO, the 
deuteron quadrupole form factor up to N 2 LO and the isoscalar amplitude in np —>■ dj up to 
N 2 LO. 

The lagrange density describing such interaction is 0] 
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We have not shown higher dimension operators, such as those corresponding to 0(p 8 ), but 
it is obvious how to include them and what the notation is. The tree-level amplitude for an 
3 Si — > 3 D\ transition resulting from this Lagrange density is (in momentum space) 
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The coefficients that appear in eq. (|T9| ) themselves have an expansion in powers of Q, e.g. 
j(j( sd ) = ic^ d \ + ^Cqq 1 + 

In order to calculate the 3 Si — 3 Di mixing parameter E\ up to N 3 LO we only require one 
insertion of dressed with the appropriate 3 Si — 3 S\ interactions. The expression for e\ 

is straightforward but long and so we do not present it here. We define the coefficients that 



arise in the momentum expansion of e±, E± and E^' , by 
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The coefficients are fit to the Nijmegen partial wave analysis []43 |, and are found to be 



E[ = 0.386 fm 2 and E\ = —2.800 fm 4 . As far as the power counting is concerned it 
is important to note that the coefficients e[ 2 ^ and E^f 1 are set by physics at the high scale 
and therefore are Q° or higher. Thus, contributions of order Q~ l must identically vanish 
H| . The superscript on the E± 's denotes the lowest order in Q at which contributions 
may arise. These conditions ensure non-trivial relations between the coefficients in eq. ([19]) 
as the renormalization scale is reduced below the high scale. Each of the coefficients in 
eq. (0) can be written in terms of physical observables, such as E[ 2 \ Ef \ p d , 7, and the 
renormalization scale [i. 

We are free to choose parameters other than the E^ 's to expand in. A quantity that 
is more directly related to the properties of the deuteron, is r] sd , (written in terms of the 
3 5i — 3 D\ mixing angle in the Blatt and Biedenharn parameterization of the S-matrix K] ) 
which is defined to be 
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evaluated at the deuteron pole, |p| = i*y. The difference between using {f] sd ,E[^} and 
{E[ 2 \ e[^} is higher order in the expansion. In terms of the coefficients E^ and e[ 4 ^ it is 
easy to show that, up to N 3 LO 
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which is, order by order, r] sd = 0.0207 - 0.0042 + 0.0076 - 0.0017 + ... . Numerically, it 

is clear that the expansion is converging, but slowly due to the relatively large size of e[^ 

(2) 

compared to E\ . This slow convergence will give rise to slow convergence in observables 
involving the deuteron and therefore it is convenient to invert this relation and use the very 
precise [|^] determination of rj sd = 0.02543 ± 0.00007 as one of the expansion parameters. 
At LO in EFT(-fr) the deuteron quadrupole moment is given entirely in terms[| of 7] sd 
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At NLO the four-nucleon-one-photon operator 
grange density 



with coefficient tJ Q defined by the La- 
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where E is the electric field operator, contributes to the deuteron electric quadrupole moment 
as well as to the isoscalar E2 amplitude in np — > dj. The counterterm is determined by 



1 In [Q] we wrote the quadrupole moment in terms of the 's. This leads to slight numerical 
differences between the two expressions (formally higher order differences). 
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fitting the NLO amplitude to the observed quadrupole moment, and it is convenient to 
define the quantity 6/j,q 

8 f i Q = ti f Q ,t '-VQ° ) = -0.0492 fm 2 , (27) 

which is taken to scale as Q 1 in the power counting. Solving for the quadrupole counterterm, 
one finds 

° Q = 71 2 7 2 (/,-7) 2 ' (8) 

Naively, higher order quadrupole counterterms contribute to the quadrupole form factor and 
quadrupole moment at N 2 LO. However, an RG analysis of such contribution shows that they 
first contribute at N 3 LO, and we can neglect them in our analysis. Explicit calculation of the 
deuteron quadrupole form factor pl|J^] up to N 2 LO gives (neglecting relativistic corrections 




that are suppressed by additional factors of m 2 /M^) 

ik Fam = ^ " [ 4|k| ( 7 (1 + + ^ + + ~<"" ] ~ \~< rW ) 

- (3|k| 2 + 16 7 2 ) ((l + IPi + l 2 pfj - tan" 

and the deuteron quadrupole moment, \iq = Fg(0)/Mj, is reproduced straightforwardly, 
r^v = 0.79 fm is the isoscalar nucleon charge radius that first enters at N 2 LO. It is clear 
that up to N 2 LO the quadrupole form factor has a well behaved expansion in powers of •ypd- 
The overall normalization is largely determined by 77,^, with the counterterm appearing at 
NLO required to reproduce the quadrupole moment @]. It is interesting to note that even 
at N 2 LO there is no contribution from and the form factor is given entirely in terms 

of r] s d 7, Pd and 5pQ. A plot of the quadrupole form factor at LO, NLO and N 2 LO can be 
found in fig. ([[]). 

Using the above analysis we are in a position to make a prediction for the isoscalar 
amplitude in np — ► dj up to N 2 LO. The LO, NLO and N 2 LO contributions to X^2 S are 



X^fJ = ~-Vsd = -2.54xl0- a 
X E2, = — gg- Vsd + — y= 5(i Q = -0.86 x 10 



1 



Z£s m = ^ + % E[ 4) = -0.67 x 10- 3 . (30) 



64 Vsd + 40 

It is clear that the perturbative expansion is converging, however, the ratio of the third to 
second term in the expansion is not particularly small. This suggests that a N 3 LO calculation 
is required before one has confidence in the value of this amplitude. A conservative estimate 
of the uncertainty in the N 2 LO calculation is the size of the N 2 LO contribution itself. 

Numerically, the EFT (ft) calculations of the subleading amplitudes for near threshold 
np — > g?7 capture areQ 



2 Our definition of the E2$ amplitude is of opposite sign to that used in [fiCfl , and hence 

7^E2 = — Xe2 s /Xm1 v 



9 



100 200 300 400 

k (MeV) 

FIG. 1. The deuteron quadrupole form factor. The dashed, dotted and solid curves denote the 
LO, NLO and PPLO deuteron quadrupole form factors computed in EFT^). The dot-dashed curve 
corresponds to a calculation with the Bonn-B potential in the formulation of JffiJ . 
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Xmi v Xmi v 

with uncertainties that we naively estimate to be of order ~ 60% and ~ 15% respectively, 
due to the omission of higher order terms. For an incident neutron speed of |v| m/s in the 
proton rest frame, we find 



Xei v 

X M i v ~ ' 1 2200 



1.2 x 10- 4 ( -LLJ- ) , (32) 



with an uncertainty that we estimate to be of order ~ 3% ||41|| . Even for neutrons with 



|v| = 2200 m/s the Ely capture cross section is comparable to the suppressed amplitudes 
for Mis an d E2 S capture. 

Using these amplitudes to compute the photon polarizations P 7 , we find 

P 7 (M1) = -7.1 x 10" 4 , P 7 (P2) = -3.5 x 10" 4 , (33) 

giving a total of P 7 = —1.06 x 10~ 3 in the forward direction, approximately 2/3 of the 
experimentally determined value of |7J] P^ x p* = —(1.5 ± 0.3) x 10 -3 . Given the large 
uncertainty in the calculation of the Mlg amplitude, and the uncertainty of the mea- 
surement, the two are consistent at the order to which we have calculated. Our value 
of P 7 (M1) = —7.1 x 10~ 4 is in complete agreement with the results of Burichenko and 
Kriplovich || of P 7 (M1) = —7.0 x 10~ 4 from a Reid soft-core calculation, but is some- 
what less than their zero-range calculation of P 7 (M1) = —9.2 x 10~ 4 . However, given 
the large uncertainty in our Mis amplitude, both values are consistent. Our value of 
P 7 (P2) = —3.5 x 10~ 4 agrees wellR with the recent calculation of Park, Kubodera, Min 



3 Given that our NLO calculation reproduces the numerical value of both the Ml$ and E2$ 
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and Rho P0| , and lies somewhere between the zero-range approximation calculation of 
P 1 (E2) = —2.4 x (which we reproduce at LO in EFT(tt)) and Reid soft-core calculation 
of Pj(E2) = -3.7 x 1(T 4 by Burichenko and Kriplovich [g. 

The power of effective field theory is that there are well-defined expansion parameters, 
even when loop graphs appear. It is therefore natural to understand the power counting 
of the spin-dependent asymmetries that we have considered. The Mis amplitude starts at 
order Q°, but receives its first non-zero contribution at order Q 1 . We have only computed 
the order Q 1 contribution. In contrast, the E2g amplitude starts at order Q 2 and we 
have computed the order Q 2 , Q 3 and Q A contributions. Therefore, the observable P 7 has 
been computed only to order Q 1 , despite our calculation of part of the order Q 2 , Q 3 and 
Q 4 contributions from the E2$ amplitude. This is apparent in the size of the uncertainty 
arising from higher order terms in the Mis amplitude, that we have discussed extensively. A 
similar statement can be made about the angular asymmetry, in particular S^ 2 \ which starts 
at order Q 2 with the interference between MI5 and E2s starting at Q 3 . Experimentally, 
measurement of both asymmetries will allow for an extraction of both Mis and E2s (when 
Ely is negligible and noting that the amplitudes are real at threshold), as is clear from 
eq. (D and eq. (§). 

In conclusion, we have used the effective field theory without pions that describes the 
nucleon-nucleon interaction to find analytic expressions for the isoscalar Ml and isoscalar 
E2 contributions to the np —>■ d^j capture process near zero incident nucleon momentum. 
The E2s amplitude is determined at the 15% level, and we find a value consistent with 
previous calculations. Due to the vanishing contribution of the one-body operator up to 
NLO, the uncertainty in the Mis amplitude is estimated to be at the 60% level. This 
relatively large uncertainty at NLO is consistent with the range of amplitudes determined 
with other approaches. A N 2 LO calculation may be able to reduce this uncertainty. However, 
additional counterterms that may arise at N 2 LO must be determined elsewhere, otherwise 
more precise predictions for these subleading amplitudes will not be possible. 



We would like to thank Mannque Rho for challenging us to make estimates of these 
suppressed amplitudes. We would also like to thank Jim Friar for a question he asked us 
two years ago about the prediction for 7] S( i in effective field theory. This work is supported 
in part by the U.S. Dept. of Energy under Grants No. DE-FG03-97ER4014. 



amplitudes computed by Park, Kubodera, Min and Rho [J10[|, the "Rho-Challenge" has been met. 
These observables do not distinguish between the two approaches at the order to which we are 
working. 

11 



REFERENCES 



[1] H. A. Bethe, Phys. Rev. 76, 38 (1949); H. A. Bethe and C. Longmire, Phys. Rev. 77, 
647 (1950). 

[2] H. P. Noyes, Nucl. Phys. 74, 508 (1965). 

[3] A.E. Cox, S.A.R. Wynchank and C.H. Collie, Nucl. Phys. 74, 497 (1965). 
[4] J.W. Chen, G. Rupak and M.J. Savage, jtiucl-th/99020561 
[5] D. O. Riska and G.E. Brown, Phys. Lett. B 38, 193 (1972). 

[6] T.-S. Park, D.-P. Min and M. Rho, Phys. Rev. Lett. 74, 4153 (1995); Nucl. Phys. A 

646, 83 (1999). 
[7] A. N. Bazhenov et al, Phys. Lett. B 289, 17 (1992). 
[8] A. P. Burichenko and I. B. Khriplovich, Nucl. Phys. A 515, 139 (1990). 
[9] T. M. Muller, Private Communication. 
[10] S. Weinberg, Phys. Lett. B 251, 288 (1990); Nucl. Phys. B 363, 3 (1991); Phys. Lett. 

B 295, 114 (1992). 

[11] C. Ordonez and U. van Kolck, Phys. Lett. B 291, 459 (1992); C. Ordonez, L. Ray and 

U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994); Phys. Rev. C 53, 2086 (1996); U. van 

Kolck, Phys. Rev. C 49, 2932 (1994). 
[12] J. Friar, |nucl-th/9601012| ; |rmcl-th/9601013| ; Few Body Systems Suppl. 99, 1 (1996); 

|nucl-th/98040iq ; J.L. Friar, D. Huber, and U. van Kolck, Phys. Rev. C 59, 53 (1999). 
[13] T.-S. Park, D.-P. Min and M. Rho, Phys. Rev. Lett. 74, 4153 (1995) ; Nucl. Phys. A 

596, 515 (1996). 

[14] T.D. Cohen, Phys. Rev. C 55, 67 (1997). D.R. Phillips and T.D. Cohen, Phys. Lett. B 
390, 7 (1997). K.A. Scaldeferri, D.R. Phillips, C.W. Kao and T.D. Cohen, Phys. Rev. 
C 56, 679 (1997). S.R. Beane, T.D. Cohen and D.R. Phillips, |nucl-th/9709062[ ; D.R. 
Phillips, S.R. Beane and T.D. Cohen, Annals Phys. 263, 255 (1998). 

[15] M.J. Savage, Phys. Rev. C 55, 2 185 (1997). 

[16] CP. Lepage, jnucl-th/9706029], Lectures at 9th Jorge Andre Swieca Summer School: 

Particles and Fields, Sao Paulo, Brazil, Feb 1997. 
[17] M. Luke and A.V. Manohar, Phys. Rev. D 55, 4129 (1997). 
[18] D.B. Kaplan, M.J. Savage, and M.B. Wise, Nucl. Phys. B 478, 629 (1996). 
[19] D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Lett. B 424, 390 (1998); Nucl. Phys. 

B 534, 329 (1998). 
[20] U. van Kolck, Nucl. Phys. A 645, 273 (1999) 

[21] D.B. Kaplan, M.J. Savage and M.B. Wise, Phys. Rev. C 59, 617 (1999). 

[22] J.W. Chen, H. W. Griesshammer, M. J. Savage and R. P. Springer, Nucl. Phys. A 644, 

221 (1999);Nucl. Phys. A 644, 245 (1999). 
[23] J.W. Chen, ^ucT-th/ 98 10021 . 

[24] M. J. Savage and R.P Springer, Nucl. Phys. A 644, 235 (1999). 

[25] D. B. Kaplan, M. J. Savage, R. P. Springer and M. B. Wise, Phys. Lett. B 449, 1 (1999). 
[26] M. J. Savage, K. A. Scaldeferri, Mark B.Wise, |nucl-th/9811029j 

[27] T. Mehen and LW. Stewart, |nucl-th/9901064| ; |nucl-th/9809095| ; Phys. Lett. B 445, 
378 (1999). 

[28] G. Rupak and N. Shoresh, |tiucl-th/9902077| . 

[29] J. Gegelia, Phys. Lett. B 429, 227 (1998); |nucl-th/9806028| ; |nucl-th/9805008| ; 
|nucl-th/9802038l . 



12 



[30] A. K. Rajantie, Nucl. Phys. B 480, 729 (1996). 

[31] J.V. Steele and R.J. Furnstahl, Nucl. Phys. A 637, 46 (1998); Nucl. Phys. A 645, 439 
(1999). 

[32] T.D. Cohen and J.M. Hansen, Phys. Rev. C 59, 13 (1999); Emcl-th/9901065| . 

[33] T.-S. Park, K. Kubodera, D.-R Min, and M. Rho, Astrophys. Jour. 507, 443 (1998); 

Phys. Rev. C 58, R637 (1998). 
[34] X. Kong and F. Ravndal, |nucl-th/9803046| ; Phys. Lett. B 450, 320 (1999); 

[nucl-th/9902064| ; |aucl-th/9904026| . 
[35] E. Epelbaoum and U.G. Meissner, |nucl-th/9902042| ; E. Epelbaoum, W. Glockle, A. 

Kruger and Ulf-G. Meissner, Nucl. Phys. A 645, 413 (1999); E. Epelbaoum, W. Glockle 

and Ulf-G. Meissner, Phys. Lett. B 439, 1 (1998); E. Epelbaoum, W. Glockle and Ulf-G. 

Meissner, Nucl. Phys. A 637, 107 (1998). 
[36] T. Mehen, I. W. Stewart and M.B. Wise, |hep-ph/9902370 . 
[37] D. R. Phillips, S. R. Beane and M. C. Birse, |hep-ph/98i0049| . 
[38] M. N. Butler and J. W. Chen, |nucl-th/9905059| . 

[39] PF. Bedaque, H.W. Hammer and U. van Kolck, Phys. Rev. Lett. 82, 463 (1999); Phys. 
Rev. C 58, R641 (1998); PF. Bedaque and U. van Kolck, Phys. Lett. B 428, 221 (1998). 

[40] T.-S. Park, K. Kubodera, D.-P. Min and M. Rho, Talk presented by M. Rho at the 
Nuclear Physics with Effective Field Theory: 1999 workshop, INT, University of Wash- 
ington, Seattle, February 1999. frucl-th/9904053| . 

[41] J.W. Chen and M.J. Savage, |nucl-th/9907042| . 

[42] H. P. Stapp, T. J. Ypsilantis and N. Metropolis, Phys. Rev. 105, 302 (1957). 

[43] V.G.J. Stoks, R.A.M. Klomp, C.P.F. Terheggen and J.J. de Swart, Phys. Rev. C49 

(1994) 2950, |nucl-th/9406039| . 
[44] S. Fleming, T. Mehen and I. S. Stewart, |nucl-th/990605q . 
[45] J. M. Blatt and L. C. Biedenharn, Phys. Rev. 86, 399 (1952). 
[46] J. Adam Jr. , H. Goller and H. Arenhdvel, Phys. Rev. C 48, 370 (1993). 



13 



